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We consider the dynamics of a quantum phase-slip junction (QPSJ) - a dual Josephson junction - connected 
to a microwave source with frequency cUniw With respect to an ordinary Josephson junction, a QPSJ can sustain 
dual Shapiro steps, consisting of well-defined current plateaus at multiple integers of ecOmw/^r in the current- 
voltage (I-V) characteristic. The experimental observation of these plateaus has been elusive up to now. We 
argue that thermal as well as quantum fluctuations can smear the I-V characteristic considerably. In order to 
understand these effects, we study a current-biased QPSJ under microwave irradiation and connected to an 
inductive and resistive environment. We find that the effect of the fluctuations is governed by the resistance of 
the environment and by the ratio of the phase-slip energy and the inductive energy. Our results are of interest for 
experiments aiming at the observation of dual Shapiro steps in QPSJ devices for the definition of a new quantum 
current standard. 

PACS numbers: 74.50+r, 74.78Na, 85.25.-j 


I. INTRODUCTION 

The Josephson junction (JJ) is one of the most used super¬ 
conducting devices in low-temperature condensed matter ex¬ 
periments. A single JJ is the building block of various sensors 
and electronic components^"— and plays an important role in 
quantum computation and informationi^— On a more funda¬ 
mental level, JJs with small capacitance have become paradig¬ 
matic systems for studying decoherence and dissipation of a 
quantum particle coupled to the external world and for ana¬ 
lyzing the transition from quantum to classical statesjii"^^ 

Many of the JJ applications are based on the Josephson ef¬ 
fect: a Cooper-pair tunneling current Ij can flow through a 
JJ in the absence of an applied bias voltage. The amplitude 
of this supercurrent is a non-linear function of the phase dif¬ 
ference (p between the two superconductors of the junction, 
Ij = Ic sin(^). The critical current Ic is the maximum Cooper- 
pair current that can be carried by the junction. A voltage 
drop Vj = {h/2e)d<p /dt appears across the junction when (p 
changes as a function of time. The classical dynamics of 
(p is ruled by the equations of motion for a fictitious parti¬ 
cle moving in a tilted washboard potential. In particular, a 
phase-locking effect can occur when the JJ is irradiated with 
microwaves of frequency cUmw'^ Then, the so-called Shapiro 
steps of constant voltage Vj^m = tn{Ti/2e)(Omw, with m integer, 
appear in the current-voltage characteristic in addition to the 
zero-voltage supercurrent state— These steps are related 
only to the fundamental constants of physics (the Planck con¬ 
stant Ti and the electron charge e) and are used in metrology 
to define the quantum voltage standard—^— The necessary 
metrological accuracy is reached at low temperatures and us¬ 
ing junctions with large Josephson energy Ej = Oq/c/ {2n) ~ 
100 meV [<l>o = h/{2e) is the superconducting flux quan¬ 
tum] and small charging energy Ec = I2C ^ 10 neV, where 
the capacitance of the junction C plays the role of the iner¬ 
tial mass in the dynamics of the phase. Moreover, the JJ 
is typically embedded in a circuit whose resistance R < Rq, 
with Rq = hl[Ae^) = 6.45 kO the superconducting resistance 


quantum. Under these conditions, thermal and quantum fluc¬ 
tuations of the phase <p are suppressed efficiently— 

The Josephson junction has an exact dual counter¬ 
part, the so-called quantum phase-slip junction (QPSJji^t^^ 
Physical realizations of QPSJ that have been discussed 
in the literature are a single Josephson junction with 
a finite capacitance ^^^^i^^" — or a linear chain of such 
Josephson junctionS j^^i'*^*" ^ and a narrow superconducting 
nanowire— With respect to an ordinary JJ, the role of the 
phase and the charge in a QPSJ is interchanged. Specifically, 
Cooper-pair tunneling is replaced by its dual process, i.e., the 
slippage by 2n of the phase difference between two well- 
defined superconducting regions of the device. As a conse¬ 
quence, the relations governing the behavior of a QPSJ are 
exactly dual to the usual Josephson relations. The voltage 
Vj = VcSin(7rq'/e) across the QPSJ is a non-linear function 
of the charge variable q, where the critical value 14 is the 
maximum voltage that the junction can sustain. The Cooper- 
pair current Ij = dq/dt is different from zero only for time- 
dependent q. As a consequence, under microwave irradia¬ 
tion, a QPSJ should sustain a set of current steps, i.e., the dual 
Shapiro steps = '«etUmw/ 

Although the occurrence of coherent phase-slip events has 
been reported for small Josephson junctions irradiated with 
microwaves,—"— a clear experimental evidence of the dual 
steps has been elusive so far. Indeed, the dual Josephson re¬ 
lations pertain to a QPSJ with a relatively well defined charge 
q, achieved when phase-slips are produced at an appreciable 
rate, a condition which is not easily compatible with the exis¬ 
tence of a well-defined underlying superconducting state. Ac¬ 
tual realizations of a QPSJ are typically operated in a regime 
where 14 is not large, so that charge fluctuations are important, 
and may well mask the dual Shapiro steps. 

In this paper, we study the role of both thermal and quantum 
fluctuations of charge on the properties of the dual Shapiro 
steps. We present the results of a combined analytical and 
numerical analysis of a QPSJ irradiated with microwaves and 
embedded in a resistive (R) and inductive (L) electromagnetic 
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FIG. 1. (Color online) (a) Circuit of a current-biased QPSJ with 
phase-slip energy Uq = embedded in a resistive (K) and in¬ 

ductive (L) electromagnetic environment. The total bias current is the 
sum of a DC component /q and an AC one /mw(f)- (b) Dual Shapiro 
steps corresponding to /n = 0, 1, and 2 for a QPSJ in the absence of 
the environment {R = o°). Parameters: a = 1.4 and ficOmw= 1- 


environment. We will see, in particular, that an important role 
is played by the inductance L, the quantity dual to the capaci¬ 
tance C of a usual Josephson junction. By duality, we expect 
that the fluctuations of the charge q are governed by the ratio 
Uq/El of the characteristic phase-slip energy Uq = leVd {2n), 
dual to the Josephson coupling energy Ej, and the inductive 
energy El = <J>q/( 2L), dual to the charging energy Eq of a 
Josephson junction.— The larger L, the smaller El and the 
larger the ratio Uq/El, thus favoring a well-defined charge 
state of the QPSJ. Recent experiments on nanowires^^s^ and 
chains of Josephson junctions^ typically achieve Uq/El ra¬ 
tios that are of the order of 10^^ — 10^*. We will analyze 
the microwave response of a QPSJ in this regime in detail 
and study in particular the resolution and accuracy of the dual 
Shapiro steps. 


II. QUALITATIVE DISCUSSION OF THE MAIN RESULTS 

We consider a current-biased QPSJ, with critical voltage 
Vc = nUo/e, which is connected in parallel to a resistive 
{R) and inductive (L) electromagnetic environment, and is 
driven by both a DC source, /q, and an AC one /mw(0 = 
/mwCosoimwf with amplitude and microwave frequency 
CUmw [see Fig. [11^)]. This circuit is related to the voltage- 
biased version by the Thevenin-Norton theorem setting Iq = 
Vq/R and j/mwl = |Fmw|/\/R^ The results for the 

Ij-Vj curve of the QPSJ of Fig.[TJa) that will be discussed in 
the following are independent of the specific choice of the ex¬ 
ternal bias. 

Let us first consider the case when the environment is ab¬ 
sent, R -A oo, in the circuit of Fig. [11^). Then the dual 
Josephson relations describing the current-biased QPSJ can 
be straightforwardly integrated yielding 

tmwl 

yjmwj = Vc ^ Jm{a) sin [nq^/e + CUfif + moh-n^t), 

m=—CO 

( 1 ) 

where J,„ is a Bessel function of the first kind. We defined 
the parameters a — ;r/mw/(ecUmw) and cob = Ttlo/e. From 



Vj/V, 

FIG. 2. (Color online) Ij-Vj characteristic obtained from the numer¬ 
ical evaluation of Eq. l l29b in the low-conductive regime, g = 0.2. 
Parameters: hgT /Uq = 0.25, a = 1.4 and hco^v//27tUo = 1. The 
three curves are obtained using three different values of the induc¬ 
tance of the environment such that Uq/El = 0.013 (red dashed line), 
Uq/El = 0.051 (blue solid line), Uq/El = 0.141 (green dotted line). 
The inset shows the relative deviation dim = TZlj/mecOmw — 1 for the 
first step, m = 1. 


this result we see that whenever the DC bias current Iq equals 
me(i)m\n/n, the QPSJ will sustain a charge-dependent DC volt¬ 
age = VcJm{oL)sm{nqQ/e). In other words, phase-locking 
occurs, leading to the appearance of a dual Shapiro step in 
the DC current-voltage characteristic of the QPSJ which is 
located at Ij,„t = mecOmw/i^, and whose width is given by 
2VcJm{a) [see Fig. [lib)]. 

When the resistance R of the environment is restored, it is 
generally a challenging task to analyze the behavior of the 
QPSJ, due to the interplay of quantum dissipative dynam¬ 
ics, non-linearity and a time-dependent signal. The combined 
effect of the application of microwaves and the presence of 
charge fluctuations induced by the resistive-inductive envi¬ 
ronment was discussed in literature in several limits. For in¬ 
stance, assuming the charge has a classical dynamics, one can 
study a Langevin equation (see Appendix [All as discussed in 
the seminal work Ref. Quantum corrections were con¬ 
sidered in Ref. but in the absence of microwaves. A full 
quantum approach of the dynamics of the junction in the pres¬ 
ence of microwaves was discussed in Refs. [33| and [3^ using a 
density matrix approach. Thermal and quantum fluctuations 
play a crucial role as they significantly affect the shape of the 
current-voltage characteristic. 

We are now in a position to state the main results of this 
paper. We use an approach that is non-perturbative in both 
the environmental coupling strength g, defined as g = Rq/R, 
and the microwave coupling strength a, i.e., a generalization 
of the so-called P(£’)— theory^ in the dual regime As we 
will see in the following, this implies that analytical results 
can only be obtained in the limit Uq/El < 1. On the other 
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hand, this corresponds to the relevant experimental situation 
where QPSJs are studied with relatively low phase-slip rates 
and not too large inductances. In the limit Uq/El < 1, we 
find that, at the first order in Uq, the QPSJ’s current-voltage 
characteristic in the presence of microwaves can be straight¬ 
forwardly obtained from the DC result without microwaves 
(see Appendix IBT i: 


tion theory for the dissipative coupling with the external envi¬ 
ronment and for the coupling with the applied microwaves. In 
Sec. lIVl we develop the non-perturbative approach. In Sec.lVl 
we discuss the results to the leading order in the t/o expansion 
focusing on the accuracy of the dual Shapiro steps and on the 
Joule heating effects. We draw our conclusions in Sec.lVIl 


m=—oo 

in agreement with a general result recently demonstrated in 
Ref. Here is given by Eqs. (IBII) and (IB2I) and cor¬ 
responds to the voltage as obtained by means of the P{E)- 
theory, as discussed in Appendix iBl Specifically, this result 
implies that the current-voltage characteristic of a QPSJ with 
t/o /El < 1 under microwave irradiation is obtained by repli¬ 
cating the known DC characteristic of the QPSJ in the ab¬ 
sence of microwaves at the positions of the current plateaus 
Ij.m = nie(Omvil'^, which are expected for a QPSJ in the ab¬ 
sence of the environment. Equation (O was derived previ¬ 
ously in Ref. |33| for the particular case where the QPSJ is 
a Josephson junction. In this paper, we analyze the conse¬ 
quences of Eq. (|2]i in detail. In particular, we discuss the role 
of the inductance L of the external impedance in view of the 
accuracy of the Shapiro steps and consider the effect of Joule 
heating produced in the resistive environment. 

We focus on the case g < 1, for which dual Shapiro steps 
clearly appear in the I-V curve. Figure |2] displays typical 
current-voltage characteristics obtained in this situation, tak¬ 
ing g = 0.2. We see that the current-voltage characteristics are 
strongly modified in the simultaneous presence of microwaves 
and charge fluctuations induced by the environment: Rather 
than being a set of discrete steps, the QPSJ’s Ij-Vj charac¬ 
teristic is a continuous curve, connecting subsequent steps, 
bending back towards a zero-voltage state in between them. 
In other words, in the presence of microwaves, a replica of 
the relevant QPSJ’s DC current-voltage characteristic (see Ap- 
pendix Island Fig.[8ll is found for each dual Shapiro step. As 
expected, in the presence of charge fluctuations, the width of 
the steps is smaller than the value 2VcJn{oi), found for g = 0; 
also, the dual steps are no longer strictly horizontal but ac¬ 
quire a small but finite linear slope. Note the role played by 
the inductance L, which limits the effects of the charge fluctu¬ 
ations. As it is clearly seen in Fig. |2] the larger L, the larger 
the width of the steps and the smaller their slopes. This can 
be seen in particular in the inset of Fig. |2] which presents the 
relative accuracy 8lm = — = Klj/mecOmvj— 1 for 

the first Shapiro step, m = 1. The inset also shows that the 
accuracy of the dual step is not only limited by charge fluctu¬ 
ations but also by a systematic shift of the step position, down 
by about 0.0015 in relative accuracy. This is due to the finite 
overlap of the various replicas. The shift can be reduced by 
increasing the microwave frequency so that the replicas are 
more separated along the /y-axis, thus reducing their overlap. 

The rest of the paper is structured as follows. In Sec. [Bll 
we introduce the model Hamiltonian for a QPSJ connected to 
a microwave source. We also show the results of the perturba- 


III. CURRENT-BIASED QPSJ 
A. QPSJ Hamiltonian 


The Hamiltonian of the current-biased QPSJ in the circuit 
depicted in Fig.[TJa) is given by 


H = —Uqcos 


■ {q + Qrl) 


m) 

2e 


^ +He, 


[{Qi},{n}] ■ 

(3) 


Here the charge and phase operators q and (p are canonically 
conjugate, satisfying the commutation relation [(p,q] = 2ei. 
As a consequence, q satisfies the equation of motion q = lit) 
and thus corresponds to the total charge injected into the par¬ 
allel combination of the QPSJ and the R-L environment. The 
first term in Eq. Q describes the nonlinear QPSJ with phase- 
slip energy t/o, which carries the charge q + Qrl, where the 
charge variable Qrl = La Qx accounts for the charge of the 
dissipative R-L environment. We model it using an infinite en¬ 
semble of harmonic oscillators (Caldeira-Leggett model) /^i^^ 
described by the third term of Hamiltonian ([3]), 


f^env [{Qa},{^a}] 


2C;i 2L;i \2e ) 


(4) 


The charge Qx and the phase ^x represent the momen¬ 
tum and position, respectively, of the A-oscillator with char¬ 
acteristic frequency (Ox = I jQLxCx- According to the 
fluctuation-dissipation theorem, (l/2)([/RL(f), fRL(0)]+) (0 — 
hco3ie[Y {(o]coth{Ticop/2), where /rl = Qrl is the fluctuat¬ 
ing current in the R-L environment and denotes the 

anticommutator. This yields the relation 


9Ie[y(®)] = 

X 





(5) 


linking the parameters of the Caldeira-Leggett bath with the 
environmental admittance of the circuit in Fig. [Ha): 


Y{(o)=l/{R-i(oL). (6) 


Finally, the coupling between the charge operator q and the 
bias current I{t) is given by the second term in Q. 

Hamiltonian Q has been used to describe QPSJs based 
on nanowires^^ Josephson junctions^ and chains of Joseph¬ 
son junctions.— In Appendix |C] we show how Eq. Q can 
be obtained starting from the well-known Hamiltonian of a 
current-biased single Josephson junction connected to a R-L 
impedance. 
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B. Current-Voltage Characteristic 

The DC current Ij flowing through the QPSJ element is 
given by the difference between the total DC current /q and 
the current flowing through the R-L impedance of the circuit 
ofFig.Ea), 

Ij = Io-Vj/R. (7) 



Here, Vj is the DC component of the voltage drop across the 
QPSJ element. Using the Josephson relation between ^ and 
Vj and the Heisenberg equation of motion for the operator 0 
generated by the Hamiltonian H, this potential reads as 


Vj 


Ti 

Ye 



= Vr(sin 


L e 


{q + Qrl) 



(8) 


The symbol (...) denotes the quantum statistical average for 
the system described by the Hamiltonian H, Eq.©. 


1. Dual Shapiro steps in the absence of environment 

By setting Qrl = 0 in Eq. Q, the coupling with the en¬ 
vironment vanishes and the system corresponds to an ideal 
current-biased QPSJ whose Hamiltonian Hq contains only the 
first two terms of ft. Introducing a complete set of discrete 
phase-states for the QPSJ, |0„) = 2n\n) with n integer, we 
can express Hq as 

Ho = --^^(|«)(u + l|+H.c.)-^^^27rn|n)(n| , (9) 

in the phase representation. When /^w = 0, Eq. (|9l) corre¬ 
sponds to the well-known Wannier-Stark ladder problem for 
a particle moving in a tilted tight-binding lattice (see Eig.[3). 
The tilt /q provides an energy difference equal to Tiojb between 
two adjacent phase states. The term proportional to Uq in¬ 
duces transitions between adjacent phase states, i.e., phase- 
slip events. In the absence of microwaves or a coupling to the 
environment, we have only coherent Bloch oscillations and 
the associated energy difference Hcob can not be accommo¬ 
dated by the system. Hence no finite DC component is found 
for the voltage Vj in this case. 

Switching on the microwave field, the tilted lattice acquires 
an additional, oscillatory slope with amplitude /mw ^ 0. Eor 
this problem, the unitary evolution operator can be evaluated 
exactly and it reads as^i^ 

U{t) = ^ (iQ) 

in which we set 


^(f) = 0)^1 +«sin (atmwO • (11) 

In Eq. (dOl, we also introduced the number operator h = 
Y,„ n \n) {n\ and the ladder operator K = Ln I”) ('^ + 11- After 
some algebra, the expectation value of the voltage operator in 


FIG. 3. (Color online) Wannier-Stark ladder. The tilt provided by the 
bias current Iq induces an energy separation hcog between adjacent 
phase states indicated by red horizontal bars, (a) Phase-locking oc¬ 
curs when the resonant condition tOg = nico^fj is satisfied. For m = 1, 
a photon with energy hcOmv/ is exchanged with the microwave source, 
(b) Environment-assisted transitions between adjacent states in the 
Wannier-Stark ladder lead also to the appearance of a finite voltage 
across the QPSJ element, (c) Wannier-Stark ladder in the presence 
of both microwave and environmental photons with energies lifBmw 
and AE respectively. 


Eq. ® on the state U{t) |^o). the time evolved initial quasi¬ 
charge state Iq'o), is 

vjmw) ^ [TTqo/e + ^(t)] 

= Vc 52 -^m (a) sin + CObI + mCOmwt) ■ 

m=-co s e / 

( 12 ) 

Equation (fTST i coincides with Eq. ([T]i and describes the ideal 
dual Shapiro steps: a non-vanishing DC-voltage now ap¬ 
pears each time the bias-current Iq = Ij satisfies the condition 
Ij = mecOmv/l as shown in Eig.dtb). The dual Shapiro steps, 
labeled with the index m = 0, ± 1,..., are replicas of the zero- 
voltage state obtained with m = 0 and a = 0, rescaled with the 
corresponding Bessel function of the first kind (a). The co¬ 
herent emission/absorption of microwave photons with energy 
hcOmv/ is at the origin of this phenomenon. The local phase 
states undergo a coherent quantum tunneling upon exchang¬ 
ing the energy JicOmw with the microwave field [see Eig.[2a)]. 


2. Perturbation theory 

We next analyze the current-voltage characteristic of the 
QPSJ in terms of perturbation theory in microwave interac¬ 
tion a and dissipative coupling g. We show that this approach 
systematically leads to divergent behavior. Eor simplicity, we 
assume the bath to be at zero temperature. 

Applying the unitary transformation = 

exp [—i0QRL/'2e] to Hamiltonian Q, we obtain the 
QPSJ Hamiltonian in the form h[ = Hq-\- Him in which we 
consider as the unperturbed Hamiltonian 

Ho = -Uocos(J^q^-^(p, (13) 

and the interaction term 

zv hi zv zv 

fiint — ^env [{Ga }; “1“ 0}] ■ (1^) 
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In this canonical form, the voltage across the QPSJ is given 
by 


= (15) 

Using the interaction picture, we expand the unitary time 
evolution operator in terms of Hint [Eq. dl] to calculate 
Vj [Eq. mi After some algebra, for vanishing microwave 
strength a = 0, we obtain for the DC component of the volt¬ 
age 


=gV^/{2RQh). (16) 

This result is linear in g and corresponds to the first order ex¬ 
pansion of the solution of the classical Langevin equation for 
the charge [see Eq. (IA3b in Appendix [Al. Its validity requires 
E//Vc ^ E hence /q ^ gVc/RQ- We conclude that perturba¬ 
tion theory breaks down in the limit of vanishing DC current 
bias. 

In the presence of microwaves, a 0, the result ([Thl l gen¬ 
eralizes to 


^(mw) _ g^c ^ /jy-, 

/o + me(Omw /’ 

which shows that the divergent behavior found for /q —?> 0 is 
repeated at the positions /q at which the dual 

Shapiro steps are expected. 

Although the perturbative approach is divergent and is in¬ 
appropriate to describe the dual Shapiro steps in the pres¬ 
ence of dissipation, it is useful for giving a simple picture of 
the QPSJ’s dynamics: The incoherent tunneling of the local¬ 
ized phase states in the Wannier-Stark ladder generally oc¬ 
curs via the combined emission and/or absorption of a cer¬ 
tain number of photons with energy /iWmw of the microwave 
source and the exchange of an amount of energy AE with 
the thermal bath (see Eig. |2l. One expects that the interplay 
between the photon-assisted and environment-assisted phase- 
slippage causes the smearing of the ideal dual Shapiro steps. 
Indeed, the sharp resonance condition Tkob = associ¬ 

ated to the single microwave photon emission can not be ful¬ 
filled anymore as the QPS junction can now dissipate the en¬ 
ergy Ticob at any bias current because the energy difference 
AE = h{(OB — Wmw) is emitted in the environment (see Pig. ID. 

IV. DUAL SHAPIRO’S STEPS: NON-PERTURBATIVE 
APPROACH 


We now develop a theory to describe the combined effect of 
charge fluctuations induced by the environment on one hand 
and microwave irradiation on the other hand, which is non- 
perturbative in both g and a. To determine the current-voltage 
characteristic of the QPSJ by means of Eq. 0, we need the 
DC component Vj of the potential across the QPSJ given by 
Eq. 0. The quantum statistical average in the right-hand 
side of Eq. 0 can be calculated by means of the Keldysh 
formalismi^i^“— Introducing the so-called Keldysh closed¬ 
time contour which goes from t = —°o to to and back to 


t = and treating the cosine term in Eq. 0 as a pertur¬ 
bation, one can obtain an exact series expansion in terms of 
the QPSJ energy I/q. In this expansion, the contribution of 
the oscillators forming the harmonic bath decouples from the 
contribution of the QPSJ charge q so that one can evaluate the 
quantum statistical averages exactly to each order. We have 
generalized this solution taking into account the presence of 
the microwave signal. The time-dependent voltage across the 
QPSJ reads as 


Vj{to) ^{-\Y(Uo\ 


Vc 2i 

^ n=0 


2n-\-l 


£ /*!.../ *2„+lEenvEq, 


where the term F^nv, 

Fenv = 

k=l 


(18) 


jfe-i 

Tlk'Mfk'-tk) 

k'=Q 


(19) 

accounts for the environment-assisted phase-slip events and 


Eq = = exp |/ £ [coBrjktk + «sin {cOmvtrjktk)] | 

( 20 ) 

is related only to the free dynamics of the charge q as given 
by Eq. ([TTIi . The dichotomic variables rjh = ±1, with k = 
0,1,... ,2n + 1, satisfy the constraint fl/t = 0 and the 

sum over all the possible configurations of 7]^. stands 

for the product of the 2n -|- 2 sums Lt]i=± ''' L 7 j,„+i=±- 

The functions of time M{t) and A(f) in Eq. ( [T^ describe 
the exchange of energy between the QPSJ and the external 
electromagnetic environment. They determine 

y(f) = —M(f) — i sign(f)A(f), (21) 

i.e., the charge-charge correlation function 

J{t) = Y{Qx{t)Qx{0)-&iiY)) , (22) 

X 

which quantifies the fluctuations of the tunneling phase due 
to the thermal bath .^’^'^^ In particular, J{t) gives the coupling 
strength between the QPSJ and the environment. Eor the 
current-biased configuration of Eig.dJa), we have 

dco ( — 1 \ 

= """e L ’ (23) 

with j3 = 1 /ksT the inverse temperature. An exact calculation 
yields^ 


A(f) = 7 rg( 


1 _ 


(24) 


M{t) =2g—— —A(f)cot —-— ) +2g Y - -^ 

np V 2 ; 1-(v„/®rl)2 

( 25 ) 
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Here v„ = 2nn/hfi is the n-th Matsubara frequency, and 
CUrl = R/Lis the frequency scale characterizing the environ¬ 
ment fluctuations at vanishing temperature. 

The Jacobi-Anger expansion exp[/a sin(x)] = 
T,m=-ooJm{ci)s^p[imx] allows to cast Fq in terms of the 
Bessel functions of the first kind Jm{<x), 




+ +0O 

^Jmoicc). ■ ■ ^ ■4i12„+i 


mo=-«> m2„+i =-~ 


(a)exp 


2h+1 

k=0 


(26) 

Performing the change of variables = h-\— 4, each time 4 
can be expressed as 4 = 4 — La=i "4 k > 1. Then Eq. (|2^ 
becomes 


2h+1 


= L '®niw4 Tlkmk 

{m^} \ mk J \ k={) j 


X exp 


2 h+1 k 

^ T (khnv/kklk') Y^ 

k=0 h=\ 


(27) 


where we used the sum rule r\k = 0. Unlike the functions 
M( 4 / — 4 ) and A( 4 / — 4 ) in Eq. ( fT9] l which depend only on 
the time difference — 4 = Y!h= \ '^h — L/,'=i Efl- <1^ is a 
function of the time 4 at which we calculate the voltage across 
the QPSJ. Prom Eq. dZTl) we observe that the frequency spec¬ 
trum of Eq. ( fTSl l at the time 4 involves integer components 
of the single fundamental frequency Omw applied to the dual 
junction. This frequency mixing is due to the QPSJ which 
operates as a non-linear capacitance, i.e., it is related to the 
cosine dependence of the QPSJ energy as a function of the 
charge q. Thus, in the steady state regime, we can extract the 
DC component by considering the time average of the general 
signal as /(f) = (1 /Tjnw) dt f{t) over a microwave pe¬ 

riod Tjnw = '27il(i)m\n where f, is an arbitrary initial time. Then, 
the DC voltage reads as 


Vj V(4) 


.••Eq(4) 





Tlk'^k ^ 


(28) 

The latter quantity is different from zero only if the sum rule 
^?=o* ^kkrik = 0 is satisfied for each arbitrary configuration of 
the variables {rjk} at given set of the integers { 111 ^} associated 
to the expansion of the Bessel functions. 


V. LOWEST ORDER RESULTS 


A general analysis of the t/o-expansion Eq. ( l28b is only pos¬ 
sible in limiting cases. We focus here on the experimentally 
most relevant limit of relatively small QPSJ energy Uq, typi¬ 
cally encountered in Josephson junction-based QPSJs. Then 
Eq. ( |28] | can be approximated with its first term. We discuss 
the range of validity of this approximation in the following. 


A. Microwave irradiated QPSJ 

Considering n = 0 only, the non-zero dichotomic variables 
are T 7 o = i and rji = ±. Since they have to satisfy the con¬ 
straint Y,k fl/t = flo + fli = it follows that the allowed config¬ 
urations {rjk} = (tjo, J?i) are and (-f,-), i.e., 770 and 771 

have opposite sign. This means that the time-average given by 
Eq. (|28 T i is different from zero if the indices niQ and mi of the 
two possible sums of Bessel functions in Eq. (l26l l are equal. 
Then Eq. (|28]) can be written as 

V^^'^\cOb)= £ 72 (a) Ej°^) (075, (29) 

m=—oo 

where voltage across the QPSJ in the absence 

of the microwaves [see Eq. (IBII) of Appendix |B]. Thus, un¬ 
der the effect of the microwave radiation, the first-order volt¬ 
age across the QPS junction is the superposition of an infi¬ 
nite number of zero-microwave potentials shifted by an inte¬ 
ger multiple m of 07mw Unlike Eq. ([T2l i. the weight of the m-th 
term in Eq. (|2^ is determined by the squared first-kind Bessel 
function of the m-th order, J^{(x). This result is in agreement 
with the general theorem proved in Ref. |5^ and was previously 
reported in Ref.lUfor the particular case of a Josephson junc¬ 
tion. Since the sum rule (a) = 1 holds, the larger is 

a, the smaller is the amplitude of the voltage corresponding 
to m = 0 and consequently the more important is the contri¬ 
bution of the higher-order terms. In other words, changing the 
amplitude a, the constant total weight re-distributes among 
the infinite terms of Eq. ( [29] ). 

Using Eq. (|7]) in combination with Eq. ( |29] |, we And that the 
Ij-Vj characteristic of the QPSJ consists of (7i7a7mw)-shifted 
and rescaled copies of the QPS J’s characteristic in the absence 
of microwaves [Eq. (IB 11) 1 obtained for /mw = 0. These fea¬ 
tures occurring at /y_„, = ineCOmwI'k^ represent the dual Shapiro 
steps smeared by quantum and thermal fluctuations induced 
by the thermal bath. These results are shown in Eig. |2l ob¬ 
tained by direct numerical evaluation of Eq. dB in com¬ 
bination with Eq. ( i29b for g < 1. The plotted smeared Ij- 
Vj curves result from the competition and interference be¬ 
tween the environment-assisted phase slippage and the pure 
photon-assisted tunneling of the phase induced by the mi¬ 
crowave held. In order for these features to be resolved, 
the microwave frequency Omw has to be much larger than 
0 )™“ Ki2ngl{Tifi ), the bias current corresponding to the back- 

bending point {vj^^/Vc)max [see Eq. (IB51 l1. 

When g > 1, the current-voltage characteristics of the 
microwave-irradiated QPSJ typically look like the ones plot¬ 
ted in Pig. m We And that they consist of replicas of the 
smeared current-voltage characteristics for g > 1 and Imw = 0 
(see Pig. [B centered around the positions of the ideal dual 
Shapiro steps. Since the Ij-Vj characteristics for g > 1 are 
more smeared than the ones found in the low-conductive 
case, a higher microwave frequency hatmxn = 20 has 
been used to resolve the various replicas in Pig. a When 
increasing the inductance L for g > 1, the smearing effects 
are reduced. The inset of Pig. a shows the relative accuracy 
5lm = 7tIj/me(Omw — 1 of the structure found at m = 1 when 
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FIG. 4. (Color online) Ij-Vj characteristics obtained from the numer¬ 
ical evaluation of Eq. l l29t in the high-conductive regime, g = 4.2. 
Parameters: kgT/U q = 0.25, a = 1.4, and hcOmv/I'^^U q = 20. The 
three curves are obtained using three different values of the induc¬ 
tance such that Uq/Ei = 0.013 (red dashed line), Uq/Ei = 0.051 
(blue solid line), Uq /£/, = 0.141 (green dotted line). The inset shows 
the relative deviation 8lm = kIj/ rnecOmvi — 1 of the structure found 
for m = 1 with respect to a perfect, first step obtained for m = 1. 


compared to a perfect, dual step. We see that the high con¬ 
ductance case does not produce single dual steps, but rather a 
doublet of two steps, located symmetrically around the value 
meaimw/?!- Combining Eq. ( |29] | and the asymptotic result 
(IB6b . we expect the positions of the steps of the doublets to 
approach their asymptotic values meoimw/± 4>o/2L with in¬ 
creasing conductance g. Eventually, a single dual step is re¬ 
covered for L oo. 


B. Accuracy of the current steps 

The reduction of quantum and thermal fluctuations affect¬ 
ing the current steps is crucial for their experimental observa¬ 
tion as well as their potential applications, such as in metrol¬ 
ogy. In this respect, it is important to analyze the accuracy 
of the dual steps. We focus on the relevant regime of low 
conductance, g < 1, where actual well-defined dual Shapiro 
steps are found and examine the smearing of the m-th step 
by considering the relative deviation 5/„, = {Ij — Ij,m)llj,m = 
nij/me(Omv/ — 1. Based on the asymptotic results of Eqs. (IB4b 
and (IB5b . we expect a minimal smearing when T and g are 
chosen as small as possible and L large. 

The behavior of 5lm as a function of some of the rele¬ 
vant system parameters is studied numerically in Figs. |5]and 
|6]for the first dual Shapiro step, m= 1. In these figures, the 
solid, dashed and dotted lines correspond to three different mi¬ 
crowave strengths a = 1.4, 2.2, and 3.2. Also shown (dashed- 
dotted line) is the behavior of the unperturbed dual Shapiro 
step for a = 2.2, i.e., 7^(2.2) x vj^^\c0B — CUmw), obtained 




FIG. 5. (Color online) Relative deviation Sim = Ttlj/mecOmw — 1 for 
the first current step, m = 1, for kgT/U q = 0.1 and /iTlUo = 
0.16. The ratio Uq/Ei is equal to 0.0032 for the plotted curves in 
panel (a) and 0.013 for the ones in (b). In both panels, the I-V char¬ 
acteristics have been obtained using three values of the microwave 
strength a: 1.4 (red dashed line), 2.2 (blue solid and dotted-dashed 
lines) and 3.2 (green dotted line). The (blue) dotted-dashed line cor¬ 
responds to the unperturbed current step for a = 2.2 (see text). 


by subtracting the contributions from all the other steps corre¬ 
sponding to m 7 ^ 1 from the signal. 

One sees that two phenomena generally limit the accuracy 
of the steps: (i) they are smeared around the actual plateau 
value and (ii) their position is offset with respect to the ex¬ 
pected one. The latter phenomenon is absent for the unper¬ 
turbed step: indeed the shift of the step position is due to the 
finite overlap of the m = 1 replica of the Bloch nose with all 
the other replicas niy^ 1. This suggests that increasing the mi¬ 
crowave frequency should yield a better accuracy of the step 
position as it separates the replicas more, thereby reducing 
their overlap and, at the same time, improving their individ- 
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FIG. 6. (Color online) Relative deviation Sl,„ = Ttlj/mecOmw — 1 for 
the first Shapiro step, m=l, forksT /(/q = 0.1, /2tzUo = 2 and 

Uq/Ei^ = 0.013. The plotted I-V characteristics have been obtained 
using three values of the microwave strength a: 1.4 (red dashed 
line), 2.2 (blue solid and dotted-dashed lines) and 3.2 (green dotted 
line). The (blue) dotted-dashed line corresponds to the unperturbed 
Shapiro step for a = 2.2 (see text). The inset shows a close view of 
the steps plotted in the main panel. 


ual resolution. The result of an increasing of (dmw on the step 
position can be seen by comparing Fig. |5lb) with Fig. |6] We 
notice, for instance, that when a = 2.2 the relative offset re¬ 
duced from about 0.02 in the former to about 0.0004 in the 
latter by increasing Omw by a factor of 10. 

It is interesting to investigate why the curve for a = 2.2 
is less affected by the offset than the one for a = 1.4, al¬ 
though the step size is the same for both curves. Indeed, the 
value of the squared Bessel functions J\{(x) determining the 
m = I step width is almost equal for the two curves. How¬ 
ever the value Tg (a) is very different; Jq{2.2) « 0.01 whereas 
7 q( 1.4) 0.32. In other words, the m = 0 dual Shapiro step 

will strongly influence the step m= I for a = 1.4, leading to 
a large offset, whereas it influences the m = I step much less 
for a = 2.2. The step corresponding to a = 3.2 is more or 
less structureless, as its weight is very small, 7^(3.2) Ri 0.07. 

As far as the smearing is concerned around the actual 
plateau position, a comparison between Figs. |3a) and |3b) 
shows the effect of the inductance. Increasing the inductance 
by a factor of 4 reduces the relative width of the step from 
about O.l in Fig. 13 a) to about 0.05 in Fig.|3b). 


C. The effect of Joule heating 

In this Section, we discuss an important aspect related to 
the experiment aimed to detect dual Shapiro steps, namely 

the effect of Joule heating in the I-V characteristic of the 
Qpcj 43.47.49.50 

As we have seen above, we expect to approach the ideal 
dual Shapiro steps of Fig. [Ttb) under the condition g <C 1. 



Vj/Vc 


(a) 



VjIVc 

(b) 

FIG. 7. (Color online) Effect of Joule heating on the dual Shapiro 
steps obtained from the numerical evaluation of Eq. l l29t in the low- 
conductive regime, g = 0.2 for (a) the current-bias circuit Eig. [TJa) 
and for (b) the equivalent voltage-hias circuit. In hoth panels, for 
the (red) dashed Ij-Vj curves the temperature is fixed to ksT/U q = 
0.25. The (blue) solid curves in (a) and (b) have been determined 
using the effective temperatures which are the solutions of 
Eqs. ([^ and 03 respectively, with fesTph/f/o = 0.25, and (/g = 
4 GHz. The electron-phonon coupling constant and the volume of 
the resistance R are Z = 10® and Q = 10^*® m^ respec- 

tively4^ All the I-V characteristics in (a) and (b) are determined set¬ 
ting U{)lEg = 0.141, flOiYaw j2711]Q = 1, and CC = 7r/mw/(ctt)n;iw) ~ 
ttVmw/(e®mw +E- coliw) = 1.4, as for the green dotted line in 
Fig. [2] The insets show the rescaled effective temperature T^g/T as 
a function of the current through the QPSJ. 


This means that the QPSJ is ideally embedded in a highly- 
dissipative environment. Such an environment is expected 
to produce also unwanted Joule heating which in turn would 
enhance the smearing of the steps. Indeed, in the low- 
conductance limit, R ^ Rq, quantum effects due to the exter¬ 
nal bath become small, whereas thermal ones induced by heat- 
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ing may become dominant. In this context, the effective elec¬ 
tronic temperature Tgff of the R-L series can be much larger 
than the phonon temperature Tph. For the circuit of Fig. [Ha), 
the current flowing through the R-L branch is Vj/R, then the 
power dissipated by the resistance is P[ = Vj /R, where Vj is 
a function of the temperature [see Eq. (l29l) 1. It follows that 
the effective temperature can be estimated by the self- 
consistent equation^ 

C0b)/{R m ■ (30) 

In this last relation, E is the material-dependent electron- 
phonon coupling constant, and D. the volume of R. Figure|71[a) 
shows the I-V curve of a QPSJ embedded in an environment 
with .g <C 1 and fixed temperature, ksT/Uo = 0.25, where 
the Joule heating is not taken into account, together with the 
dual Shapiro steps smeared by the voltage-dependent effective 
temperature Eq. (l30l l which accounts for the exchange of en¬ 
ergy between the electrons and the phonons in the resistance 
R. We notice a reduction of the width of the steps, as one 
expects. Erom the inset of Eig. |2ta), we see that T^-ff follows 
the oscillating trend of Vj. In particular, it coincides with Eph 
whenever Vj = 0 and reaches its relative maxima for the val¬ 
ues of Vj around the maximum amplitude of the steps: the 
wider are the steps in the absence of Joule heating the larger 
is their effective thermal smearing. 

On the other hand. Joule heating affects differently the dual 
Shapiro steps appearing in the I-V characteristic of a voltage- 
biased QPSJ. In this configuration, the power Py = Ij^R, 
which is dissipated by the resistance R, is determined by the 
current flowing through both R and the QPSJ, i.e., 

Ij = {Vo-Vj)/R. (31) 

Here Vj = Vj{T,Vo) is obtained from Eq. ( |29] | replacing Iq 
with Vq/R and |/mw| with As a result, 

the effective temperature of the environment can be written as 

Teff = 7pl+//(7;ff,yo)j;/(En). (32) 

Inserting into Eq. OTl i the temperatures Feff obtained by solv¬ 
ing self-consistently Eq. ( l32l i for different values of the DC 
voltage bias Vq, we obtain the (blue) solid QPSJ’s I-V char¬ 
acteristic shown in Eig. [Tjb). Notice that this curve is more 
smeared than the one found in the current-biased case and 
plotted in Eig. |7Ja) using the same set of parameters. As 
shown in the inset of Eig. |2lb). Tiff increases with |/j| and 
is equal to Fph only when Ij ~ 0. In particular, the effective 
temperature given by Eq. (l32T i is much larger than Fph when Ij 
is close to ecOm-w/ Consequently, the Joule heating affects 
the steps for m ^ 0 more than the one occurring for m = 0, 
as one can see from Eig. |2lb), thereby compromising their 
experimental observation. The reduction of this effect is pos¬ 
sible, for instance, with the decreasing of the microwave fre¬ 
quency Omw However, the use of smaller Omw leads also to 
the increasing of the offset of the steps which we discussed in 
Sec.rVBl 

In principle. Joule effect can be reduced by increasing the 
inductance L of the environment rather than the resistance R. 


L plays the same role of R in the reduction of the fluctuations, 
as shown previously. As the dual Shapiro steps are replicas 
of the I-V characteristic at low current, we can estimate the 
leading dependence for the smearing by considering Eq. (IB4I) . 
We obtain the slope GqRq k. 2g{kBT/Uo)^{Ej/Uo)^^ for g <C 
1. We observe that the smearing due to the temperature can 
partially be compensated by increasing the inductance of the 
environment. 


VI. CONCLUSIONS 

In this paper, we discussed the microwave response of a 
QPSJ embedded in an inductive-resistive environment. We 
focused on the regime of relatively small ratio of phase-slip 
energy Uo over inductive energy Ej. The response consists of 
a series of well-defined current Shapiro steps, located at multi¬ 
ples of e(Omwl if the environmental resistance is sufficiently 
large, such that the dimensionless conductance g < 1. These 
steps are in fact replicas of the QPSJ’s Bloch nose, observed 
in the absence of microwaves. Charge fluctuations induced 
by the environment smear the steps. This smearing can be 
reduced by decreasing the dimensionless environmental con¬ 
ductance g, decreasing the dimensionless temperature kBP/U q 
and increasing the ratio Uo/Ei, which can be achieved by in¬ 
creasing the environmental inductance L. Einally, we showed 
that the conductance g can not be decreased indefinitely, as 
heating effects may develop in the environment. 

The results presented in this paper are relevant for 
recent experiments on Josephson junction chains^ and 
nanowires4i^ In these works, typical phase-slip energies Uo 
are in the range of 1 — 10 GHz, whereas the environmental in¬ 
ductances L are 50 — 500 nH. This motivated the parameter 
choices used in this paper: UojEj ranges from 0.001 — 0.1; at 
typical cryostat temperatures kBP jUo ~ 0.1 — 0.2. We found 
that, although dual Shapiro-type features could be visible ex¬ 
perimentally for these parameters, their relative accuracy re¬ 
mains limited to about 0.001 by fluctuation effects. 

To date, a systematic evidence for the existence of dual 
Shapiro steps is still lacking. The reason for this might well 
be that fluctuation effects have so far masked the steps for 
QPSJs with intermediate ratios of the parameter Uq/Ej and 
not too small conductance g. Work on nanowire-based QPSJs 
with larger values of the ratio Uq/Ej and lower conductances 
g seems promising42i^ At the same time these systems suffer 
from substantial heating effects4^ We conclude that further 
work is necessary, both on nanowires and on Josephson junc¬ 
tion chains. 
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Appendix A: Classical dynamics and Langevin equation for the 
charge 

After a unitary transformation in Eq. Q, it is possible to 
show that the charge qj on the QPSJ satisfies the Fangevin 
equation 


dqj/dt ^ Iq + - lRL{t) + 5l{t), (Al) 

as discussed in previous works, e.g. in Ref. |29l Notice 
that dqj/dt = Ij corresponds to the current flowing through 
the QPSJ. The first and second term in the right-hand side 
of Eq. (lAlb give the total bias current. The third term 
is the current flowing through the resistive-inductive branch 
of the circuit /«/,(?) = Jdt'Y{t — t')Vj{t'), where Vj{t) = 
Vcsin[nqj{t)/e\ and Y{t) is the inverse Fourier transform of 
the admittance Eq. (|6ll. The last term of the right-hand side in 
Eq. (lAlb is a fluctuating current 5l{t) of zero average whose 
Fourier component satisfies the fluctuation-dissipation theo¬ 
rem 

{5l{(o)5l{(o')) = 2 k5{o} + (o')h(o^&[Y {o})\coth{rm/IksT), 

(A2) 

where T is the temperature of the environment. When /^vv = 
0, Eq. (lAlb reduces to the well-known Langevin problem of 
the quasi-charge dynamics in the overdamped regime. Then, 
the DC current-voltage characteristic of the QPSJ corresponds 
to the so-called Bloch nose (see Fig.[8]l. 

Disregarding, for instance, the fluctuation 5/ = 0 and con¬ 
sidering the limit gUo/E l 1 in Eq. (lAlb . by direct integra¬ 
tion, one obtains the following DC voltage 


y(8l=i)) 


Q ( 

g V ^ 


/ 


8 


2 


, (A3) 


where 0(U) is the Heaviside step function. The corresponding 
current-voltage characteristic is shown in the inset of Fig. [8] 
It consists of a zero-current branch at finite voltage up to Vc 
which bends back to a low-voltage, finite current branch. Fi¬ 
nite, classical charge fluctuations, 5l{t) ^ 0, prevent the for¬ 
mation of a sharp feature in the current-voltage characteristic, 
even for small g, and yielding a smearing of the Bloch nose. 


Appendix B: QPSJ in the absence of microwaves. 

In this appendix, we consider the current-voltage char¬ 
acteristics of a QPSJ without microwave irradiation using 
Eqs. (fTSI) and (l28T l. and recall some of the results provided 
by the P(£’)—theorywhich describes the phase slippage in 
the presence of an external environmenti^i^ Setting a = 0 in 
Eq. (l28b . and retaining the term n —0 only, the voltage drop 


on the QPSJ as a function of (Ob reads as 
n 

(cob) ^ xt/o lP(hC0B)-P(-rtC0B)] , (Bl) 
yc ^ 

where we defined the function^ 

1 P+oo 

P(A£) =- / dx ^ (g2) 

Inn J-oo 

The function P(AE) represents the probability density that 
the QPSJ absorbs (AE > 0) or emits {AE < 0) an amount 
of energy |A£| from or to the external environment respec¬ 
tively during a phase-slip event. We see that an incoherent 
phase slippage by = 2n in the Wannier-Stark Ladder takes 
place only if the system exchanges the energy AE = HcOb = 
(A(p)hlo/(2e) with the environment [see Fig. |3b)]. As the 
energy spectrum of the bath is continuous, the QPSJ has a dis¬ 
sipative behavior for any value of the applied DC current Iq. 



Vj IV, 


FIG. 8. (Color online) Current-voltage characteristics of a QPSJ em¬ 
bedded in an inductive-resistive environment in the absence of mi¬ 
crowaves. Parameters: fegr/f/o = 0.25, I/q/Ez, = 0.013. Curves 
from bottom to top correspond to g = 0.2, 4.2, and <». The inset 
shows the I-V curve obtained by means of Eq. l lA3b . 


We summarize the behavior of Eqs. dB and (IB2b com¬ 
bined with Eqs. (|6]l and ( l23b in the main panel of Fig. where 
we plotted the QPSJ’s current-voltage characteristic for vari¬ 
ous values of g at low temperature as obtained by direct nu¬ 
merical integration. 

The validity of Eq. (IBlb is given by the condition 
"C 1, yielding t/omax[P(/ia)B)] ^ 1.— Analytical 
results are available, for instance, in the limit of low tem¬ 
perature and small conductance so that fiEL/2K^g ^ 1 and 

j3£L/27rV » 1- Then,^ 


yjDC) 

Uc 


(®b) 


\T{g + i^tlCOB/2K)\^ 

r( 2 ^) 


sinh(j3/la)B/2) , 


(B3) 
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where u — {l5Uo/4-n){PELe'^/2n^g)^^^ with 7 = 0.577... the 
Euler constant. Hence we find a linear conductance Go at van¬ 
ishing current Ij and voltage Vj (hcOBji /27r ^ 1), given by 


GoRQ^4{Uol5f^^-^'> 


2 g 



r 2 (g)" 

(B4) 


We note that Go ~ and thus decreases with decreasing 

temperature; similarly Go ~ 1 and thus decreases 

with increasing inductance. Moreover, Go decreases with de¬ 
creasing g. Increasing (Ob until k, 2ng/{hfi), we reach 
the back-bending point corresponding to the maximum value 


^nu = {2n^gf^e-^y^ , 

(B5) 

for g <C 1. We see that the lower the temperature T, the larger 
is the inductance L, and the smaller the conductance g, the 
closer (V'j^*^^/Ve)max is to the maximum value 14. Beyond the 
back-bending point, corresponding to the maximum voltage, 
the system enters into the Bloch oscillation branch where the 
bias energy Hcob becomes dominant with respect to both quan¬ 
tum and thermal fluctuations and the DC voltage Vj decreases 
exponentially to zero. 

When the resistance R is reduced so that g > 1, the current- 
voltage characteristic is smeared into a smooth curve with a 
maximum voltage at finite current. In the high-conductance 
regime g 1, we have that the Bloch nose broadens into a 
Gaussian: 




P{AE) 


1 

y/AnEikBT 


exp{-(A£-£i)V(4£ikB7’)}. (B6) 


As a result, phase-slip events in a current-biased QPSJ can 
only occur if the energy hnlo/e exchanged with the induc¬ 
tive environment equals E^, viz., the current Ij at which the 
QPSJ sustains the largest voltage approaches the value <t>o/2E. 
This phenomenon is dual to the Coulomb blockade found 
in a Josephson junction in a highly resistive environment, 
where the voltage at which the Josephson junction sustains 
the largest current approaches the value 2e /2C.— 


Appendix C: QPSJ Hamiltonian for an underdamped 
Josephson junction 


In this appendix, we will study the Hamiltonian of the cir¬ 
cuit ofFig.|9l whichis formedby a Josephson junction (JJ),bi¬ 
ased by a time-dependent current /(f) = /q -I-4™ cos ft)mw4 in 
parallel with a capacitance C and an external electromagnetic 
environment composed by a resistance R and an inductance L 
in series. In particular, we will show that this Hamiltonian re¬ 
duces to the QPSJ Hamiltonian Q used in the main text under 
suitable conditions. 

Neglecting the contribution of the quasi-particle excita¬ 
tions, the Hamiltonian corresponding to the circuit of Fig. |9] 
is given by the sum of the charging energy of the capacitance 
C, the non-linear Josephson energy and the energy of the en¬ 
vironment. 




2 


dt'l{t') + Q^^ + Q 


Ejcos (^) Heny . 


The phase operator ^ is the phase-difference between the 
two superconductors forming the junction and Q is its con¬ 
jugate charge operator [(p,Q] = 2e!, i.e., the charge tunnel¬ 
ing through the junction. In Eq. (ICll) . we also introduced 
2rl = La Ql which accounts for the charge noise produced 
by the R-L environment, as discussed in the main text. 



FIG. 9. Current-biased Josephson junction with Josephson energy 
Ej in parallel with a capacitance C and embedded in a resistive (R) 
and inductive (L) electromagnetic environment. The circuit is biased 
with a time-dependent current I{t). 


The equivalence between Hamiltonian (IClb and QPSJ- 
Hamiltonian (O can be demonstrated through the follow¬ 
ing steps. First, we apply the gauge and the unitary 
transformations Ug{t) = exp [—i^ fdt'I{t')/2ej and Genv = 
exp [—i^QRh/'^e ], respectively, to Eq. (IClb and we get 

^ —EjCOs{0) — -^^0 + Heny + ^}] ■ 

(C2) 

Here, the first and second term correspond to the standard 
Hamiltonian Hj of an isolated JJ. In the tight-binding regime, 
Ej Ec, Hg becomes 


Hg = —Go cos 


V e 


hi{t) 

2e 


(C3) 

where q is the quasi-ch arge operator and Go = 
Sy/EjhcOp/nexp {— y/SEj/Ec) = eVc/fr the half-bandwidth 


of the first Bloch band of Hj. Within this limit, an energy 
gap of the order of the plasma frequency TicOp = y/WEjEc 
separates the first from the second Bloch band. We 
neglect the possibility of inter-band Fandau-Zener transi¬ 
tions assuming the low temperature and bias current limit 
(kBT,hlQ/2e,hIgg,yi/2e) <C licOp as well as considering an 
off-resonance microwave field, cUmw ^ ®p. 


Finally, we apply the inverse unitary transformation Ggj,(, to 
Eq. (IC3b and we obtain the effective low-energy Hamiltonian 


G=—Gqcos 




mi 

2e 


^ + He, 


(C4) 


This is the energy operator (O of the main text describ¬ 
ing a current-biased quantum phase-slip junction coupled to 
an external R-L electromagnetic environment, as depicted in 

Fig.ma). 


(Cl) 
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